Freak waves in the linear regime: A microwave study 
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Microwave transport experiments have been performed in a quasi-two-dimensional resonator with 
randomly distributed conical scatterers. At high frequencies, the flow shows branching structures 
similar to those observed in stationary imaging of electron flow. Semiclassical simulations confirm 
that caustics in the ray dynamics are responsible for these structures. At lower frequencies, large 
deviations from Rayleigh's law for the wave height distribution are observed, which can only partially 
be described by existing multiple-scattering theories. In particular there are "hot spots" with 
intensities far beyond those expected in a random wave field. The results are analogous to flow 
patterns observed in the ocean in the presence of spatially varying currents or depth variations in 
the sea floor, where branches and hot spots lead to an enhanced frequency of freak or rogue wave 
formation. 

PACS numbers: 42.25.Dd, 05.45.Mt, 42.25.Bs 



In this Letter, we present a microwave transport 
study through a scattering system composed of randomly 
placed metallic cones, each mimicking an potential 
on the scale of its radius. For wavelengths smaller than 
or comparable to the scatterer size, we find branching 
structures that are reminiscent of electron current distri- 
butions seen in two-dimensional electron systems [l[ . At 
wavelengths larger than the cone size, the bulk of the in- 
tensity distribution approaches a multiple-scattering cor- 
rection to Rayleigh statistics, as expected in multiple- 
scattering media [2]. However, the probability of find- 
ing very high intensities is still greatly enhanced. Even 
larger fluctuations are observed after Fourier transform- 
ing the frequency-domain measurements into the time- 
domain; the extreme time-domain events may be com- 
pared to freak wave events in the ocean [3] . The enhanced 
probability of large intensities at longer wavelength is a 
residual of presumably stronger enhancements at shorter 
wavelength. 

Scanning tunneling microscopy studies of electron flow 
injected through point contacts in a high- mobility two- 
dimensional electron gas (2DEG) by Topinka etal. [l| 
exhibited intricate branching patterns of fractal appear- 
ance. This behavior was in contrast to the simple random 
wave prediction for the probability distribution of wave 
intensities, PRayieigh(^) = where / ~ is normal- 
ized to one, and is the wave height. Topinka etal. 
showed that the evolution of caustics (singularities of the 
ray density) in a random potential (with a given strength 
and correlation length) is responsible for the branching 
patterns. Kaplan later showed that the probability distri- 
bution of the branches may be computed analytically Q . 
The phenomenology holds equally well for the evolution 
of wave patterns in the sea. In shallow water, wave focus- 
ing may be caused by depth variations in the sea floor, 
and may lead to amplification of tsunami waves [H, Q . In 



deep water, the effect of eddy currents has been studied 
by Heller etal. [3], who showed that even after account- 
ing for dispersion in wavelength and direction, random 
currents greatly increase the likelihood of large ampli- 
tude events; they argued that such events might act as 
a trigger for nonlinear instability effects in freak wave 
formation. Thus, without diminishing the importance 
of nonlinear processes, an understanding of the linear 
regime is essential for the proper understanding of freak 
wave physics. 




FIG. 1: (Color online) Photograph of one of the two scatter- 
ing arrangements used. The platform has width 260 mm and 
length 360 mm. Each cone has diameter 25 mm and height 
15 mm. The probe antenna is fixed in a horizontally movable 
top plate located 20mm above the bottom (not shown). 

Motivated primarily by the 2DEG experiments [1], we 
undertook a microwave experiment to study the trans- 
port of waves through an arrangement of randomly dis- 
tributed scatterers. Figure [1] shows a photograph of the 
experimental setup. The metallic bottom plate supports 
the scattering arrangements made up of about 55 to 60 
brass cones. The source antenna is mounted close to one 
of the short sides, and varying its position enables the 
incoming waves to arrive from different directions. The 
drain antenna is mounted in the top plate (not shown), 
and acts as a weak probe. The top plate can be moved in 
both horizontal directions, allowing for a spatial mapping 
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FIG. 2: (Color online) Comparison of an experimental wave 
pattern with a classical ray simulation. Left: A wavefunction 
at frequency / = 30.95 GHz. Right: The corresponding semi- 
classical simulation, with modes 1 through 4 added together. 
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of the wave fields within the scattering arrangement. 

For quasi-two-dimensional systems with parallel top 
and bottom plates separated by a vertical distance d 
(without scatterers), the electromagnetic wave equations 
reduce to a single, scalar equation for the perpendicular 
component E±{x^ z) of the electric field. Furthermore, 
we can write E±{x,y^z) = E{x^y) cos{n7rz/d), where n 
is the transverse quantum number. This results in a two- 
dimensional wave equation for E{x^y)^ 



dy^ 



\~dJ 



E{x,y) = eE{x,y), (1) 



which for n = is equivalent to the free stationary 
Schrodinger equation in the plane [7] . The number of ac- 
tive modes n depends on frequency. Eq. ([T]) remains ap- 
proximately true when d varies slowly with position (on 
the scale of the wavelength). For n > 0, the additional 
term then mimics a potential, V{x^ y) = [n7T/d{x^ y)]'^ 
Each cone in the experiment corresponds to a repul- 
sive central potential: V{r) oc [C -\- min(r, ro)]~^, where 
ro = 12.5 mm is the cone radius and C = 16.7mm is a 
constant determined by the cone geometry. At higher 
orders, however, the height variation leads to a mixing 
among the modes. 

The left panel of Fig. [2] shows a typical wave pattern 
observed for / = 30 to 40 GHz, the upper limit accessi- 
ble by our equipment. At these frequencies, modes n = 
to 4 are propagating. The wavelengths are between 7.5 
and 10mm, i.e., somewhat smaller than but comparable 
to the cone diameter. We observe exponential decay of 
the wave intensity with distance from the source, caused 
in part by escape of the waves from the scattering setup 
(since the system is open along the perimeter). This de- 
cay is suppressed in all plots. The right panel of Fig. [2] 
shows the results of a ray simulation, obtained by solving 
the classical equations of motion in the potential gener- 
ated by the scatterers. The pattern bears a striking sim- 
ilarity with the branchlike structures found by Topinka 
etal. [H in studies of 2DEG electron flow. The present 
results provide strong evidence in favor of the conjec- 
ture [H, 0] that random potentials correlated on the scale 
of a wavelength are responsible for these features. 



FIG. 3: (Color online) Probability distribution of intensities. 
The dark (black) histogram includes all data, while the light 
(yellow) histogram excludes frequencies associated with the 
hot spots. The dotted line is the Rayleigh distribution, while 
the dashed (blue) line is a best fit using the theoretical dis- 
tribution given by Eq. (O (7 23.5). 



To avoid mixing of up to five different modes, most of 
the experiments have been performed between 7.5 and 
15 GHz, where only the first two modes n = and 1 are 
propagating. In this regime, however, the wavelengths 
are large compared to the scatterer size. Hence the sep- 
aration of the z component in the wave equation is no 
longer justified, and the interpretation of the cones in 
terms of a classical particle potential becomes invalid. 

In the weakly disordered regime, Nieuwenhuizen and 
van Rossum [2] calculated the first perturbative correc- 
tion to Rayleigh's law, yielding Ppert(^) ^ e~^[l + (/^ — 
4/ + 2)/3^] for the intensity distribution, where g is the 
dimensionless conductance, exactly what had been found 
in a microwave transport study by Genack and Gar- 
cia [9] . A nonperturbative expression has been calculated 
by Mirlin et al. [10] for the transmission between two an- 
tennas embedded in a quasi-one-dimensional disordered 
sample. 



Pdis(^) - exp ■ 



In^ 




(2) 



where 7 is proportional to the conductance and depends 
also on the position of the two antennas [10]. For 7 ^ /, 
we havepdis(^) ^ e~^+^^ where the P term gives the 
first perturbative correction to Rayleigh's law. We em- 
phasize that all corrections to Rayleigh contained in such 
multiple-scattering expressions disappear in the limit of 
short wavelength (high frequency or large g). This is not 
the case with smooth potentials, which semiclassically 
deflect the flow; indeed the effect of classical caustics be- 
comes more pronounced at short wavelengths [4] . 

Figure [3] shows the intensity distribution found in our 
experiments, averaged over the complete data set (two 
scattering arrangements, three source antenna positions 
for one of the arrangements, and a frequency range of 
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FIG. 4: (Color online) A "hot spot", observed at a frequency 
of 8.85 GHz. The experimental probability density for observ- 
ing such a hot spot is one to two orders of magnitude larger 
than that expected from multiple-scattering theory. 



7.5 to 11 GHz). The distribution is well described by 
Eq. ([2j) over three orders of magnitude. Here the situa- 
tion is comparable to the one found by the Genack group 
in a number of studies 



But for the very high in- 



tensities, the probability exceeds the multiple-scattering 
theory prediction by one to two orders of magnitude. 

Examining the data, we find that just two regions (and 
another one at the border of significance) are responsible 
for these deviations, one of them shown in Fig. |4j Each 
of these "hot spots" exists only in a limited frequency 
window about 500 MHz in width, and the range of in- 
coming wave directions able to excite each hot spot is 
only about 20 degrees wide. If the frequency ranges con- 
taining these hot spots are omitted from the analysis, the 
resulting intensity distribution is in full agreement with 
the multiple-scattering prediction, see Fig. [3l 

For the study of time-dependent waves, such as those 
found in the sea, we must superimpose waves with dif- 
ferent frequencies, entering from different directions. To 
this end we concentrate on one hot spot found at 9.5 GHz 
near the center of our scattering arrangement. We vary 
the antenna position over 80 mm along the short side of 
the scattering arrangement (see Fig. [1]), and record the 
resulting field pattern for each antenna position and each 
frequency. Then time-dependent wave fields are gener- 
ated by superposition of = 150 patterns. 



N 
i=l 



(3) 



Here tpi (r ) is a wave pattern at frequency fi , excited with 
the source antenna at position x^, and (pi is a random 
phase. The randomly chosen frequencies fi are normally 
distributed with the average at 9.5 GHz and a standard 
deviation of 0.764 GHz, thus covering approximately the 
frequency window in which this hot spot is present. Sim- 
ilarly, the antenna positions Xi are taken within the angle 
of acceptance of the hot spot. 
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FIG. 5: (Color online) Intensity distribution for the time- 
dependent wave patterns generated by Eq. Q for all mea- 
sured points [dark (black) crosses], for the hot spot region 
only [light (orange) solid circles], and for all points excluding 
the hot spot [light (blue) open squares] . The dotted line is the 
random wave expectation, while the solid (red) line is given 
by Eq. JG)). The arrow indicates the extreme event studied 
more closely in Fig. [71 The inset shows a sketch of the set 
up: the dark crosses mark the different exciting antenna po- 
sitions, the (red) shaded region corresponds to the measured 
field, and the empty rectangle inside the measured field indi- 
cates the hot spot region. The cross inside the hot spot region 
is the position of the maximal measured intensity. 
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FIG. 6: (Color online) Distribution of the time-averaged in- 
tensities s found for the 780 pixels of our measurement. The 
inset shows the same data using a semilogarithmic scale. The 
solid curve is a distribution with zy = 32 degrees of freedom. 



Fixing the probe position, we always find a Rayleigh 
law for the distribution of intensities in a time sequence, 
Pioc(^) = s~^e~^/^, but with the time-averaged value 
s = (/) depending on position. This is nothing but a 
manifestation of the central limit theorem. An example 
taken at the hot spot is given by the light (orange) solid 
circles in Fig. [5l This is the situation a ship experiences 
at a given position. If now the ship changes its posi- 
tion, or, alternatively, if the random currents change, a 
Rayleigh law with another s is found. The overall distri- 
bution of time-dependent intensities, collected over po- 
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sition and/or realization of the disorder, is then given 

by 



P{I) = J dss-^g{s)e-'/r 



(4) 



where g{s) is the probabihty density to find a local time- 
averaged intensity s. 

Normalizing the intensity s to one, we find that g{s) 
can be very well described by a chi-square distribution 



1 



r(f) 



S2 



exp 



where the number of degrees of freedom u increases lin- 
early both with the frequency range used and with the 
range of source antenna positions. This is exactly what 
is expected for a pulse made up of 1^/2 independent pat- 
terns obeying Rayleigh distributions. A typical example 
is shown in Fig. [6l If the input patterns are not Rayleigh 
distributed, a correspondingly modified distribution 
is obtained. The corrections to expression ([5]) become 
significant in the range of high amplitudes. 

With expression (j5j) for g{s), the integral (j4]) yields the 
K distribution [l2j 



P(/) = 



r(|) 




(6) 



where Kjy{x) is a modified Bessel function. The solid 
(red) line in Fig. [5] is calculated from Eq. ([6]). It fits 
nicely with the intensity distribution found if the hot 
spot region is excluded (blue squares), but not with the 
distribution including the hot spot (black crosses). This 
is not really a surprise. Already the inset of Fig. [6] shows 
that the distribution, though generally working well, 
fails to describe the rare events in the high-intensity tail. 

Figure [7] shows the most extreme event found in our 
time series, marked by an arrow in Fig. [5l and the inset 
shows the entire region about the hot spot at the mo- 
ment of this freak event. In the experiment, we observe 
events of this magnitude or greater with a probability of 
1.3-10"^. Thus, such events are still quite rare, but the 
probability is enhanced by 5 orders of magnitude com- 
pared with Eq. (|6j), and by 15 orders of magnitude com- 
pared to the Rayleigh distribution! 

This work has demonstrated the virtues of microwave 
techniques for obtaining detailed information on wave 
transport through disordered surroundings. Simply by 
varying the frequency, we are able to study both ray 
dominated branching behavior of fiow in a potential 
landscape, as well as the diffr active multiple-scattering 
regime. The interpretation of the hot spots in this latter 
regime has to remain speculative for the moment, in view 
of the small number of such hot spots showing up in the 
experiments. However, the narrow angular acceptance of 
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FIG. 7: (Color online) A freak wave event. The time evolu- 
tion of wave intensity at the center of one of the hot spots is 
shown for the most extreme event observed. The inset shows 
the region surrounding the hot spot at the moment of the 
freak event. A movie of the time evolution in the entire field 
surrounding the hot spot is available [13.] . 



each hot spot, the visually obvious branching behavior at 
the higher frequencies, and particularly the fact that the 
observed deviations from Rayleigh statistics get stronger 
rather than weaker at shorter wavelengths, all support 
the hypothesis that the hot spots are not resonant wave 
phenomena but instead are remnants of singularities in 
the classical dynamics. 
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See supplementary material at |http://link.aps.org7 sup- 
plemental/10. 1103/PhysRevLett. 104.093901 for a movie 
showing a freak wave event.. 



